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We examine a scenario of the abelianized Glasma evolution with accounting for back-reaction of 
partonic medium in ultrarelativistic heavy-ion collisions. We announce that such a generalization 
leads to the instabilities and the presence of negative color conductivity in the system. 
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I. INTRODUCTION 

Phenomenological analyses of experimental data indi- 
cate that the quark-gluon plasma (QGP) can be formed 
in ultrarelativistic A+A collisions [1]. Its local thermal- 
ization and isotropization should be mainly related to the 
fast processes stimulated by instabilities at small times 
after collision 0, H| . 

In the present theoretical picture of ultrarelativistic 
heavy- ion collisions Q , the early stage is preferably char- 
acterized by the large number of partons with "small" 
momenta of the order of the so-called saturation mo- 
mentum A a , which are better viewed as a classical Yang- 
Mills field in vacuum [5[ , sometimes named as "Glasma" . 
The initial conditions for Glasma evolution are deter- 
mined by the Color Glass Condensate (CGC) concept by 
McLerran-Venugopalan (MV) @ , where the field sources 
before collision can be presented by the randomly dis- 
tributed valent quarks of colliding hadrons and are lo- 
cated (due to Lorentz contraction) on infinitesimally thin 
sheets running along the light-cone. These sources are 
also treated as the hard partons with "large" momenta, 
which escape quickly from the system after collision. 
Thus, the original MV model neglects the interaction be- 
tween the field and the hard partonic medium. 

The space-time dynamics of the Yang-Mills fields in 
vacuum ("the melting of CGC") in assumption of boost 
invariance was investigated numerically, and the energy 
and the number distributions of the classically produced 
gluons were computed (see, for example, review Q). 
Moreover, it was shown that the violations of boost in- 
variance cause a non-Abelian Weibel instability Q lead- 
ing the field (soft) modes to grow with proper time 
However, the effect of isotropization is out of this model. 

On the other hand, the hard partons (produced after 
moment of collision) with large transverse momentum pt 
can be studied within the framework of transport theory, 
and if the presence of the soft classical field is neglected, 
the time evolution of the these partons is described by 
Boltzmann equation with a collision kernel [Iol - [i~5l | (for 
comment, see [29]). However, it has been argued that 
the collective processes caused by the soft gauge field 
should be dominant in equilibration of QGP instabili- 



ties developed due to anisotropic distributions of released 
hard partons [l6l - [T9j . 

The third regime where the back-reaction of the field 
on the hard partons (treated as particles) is still weak 
but where the self-interaction of the former may be 
strongly nonlinear is governed by a "hard-loop" effective 
action which has been derived in Ref. [20| for arbitrary 
momentum-space anisotropies. 

It is interesting to note that the numerical studies 
of anisotropic hard partonic modes coupled to unsta- 
ble soft modes revealed the tendency of the non-Abelian 
gau ge fi elds to "abelianize" during the stage of instabili- 
ties [13, Hi]. It means that the field commutators become 
much smaller than the fields themselves. Moreover, the 
dynamics of the Abelian and non-Abelian fields is quali- 
tatively the same, if these fields are not strong enough. 

In this paper, we examine a scenario of Glasma evolu- 
tion with CGC-like initial conditions, when the presence 
of the (momentum-)anisotropic medium of hard partons 
is also taken into account. Our goal is to evaluate ana- 
lytically the behavior of such a system in short-time in- 
terval in weak-interaction regime, when the application 
of the abelianized version of the field dynamics is pos- 
sible. Although the last condition demands to consider 
the system at relatively large times after A+A collision 
(as follows from numerical investigations) but simplifies 
the problem considerably. However, it is already pointed 
out in Refs. f22hl2~4| that the early equilibration of QGP is 
not necessary to describe pion and kaon spectra observed 
experimentally at RHIC in Brookhaven. 

Since the momentum-space anisotropy of the system 
can be estimated by means of transport coefficients, we 
attempt to calculate a conductivity tensor and to deter- 
mine an effect of instabilities on it. It is expected that the 
back-reaction can lead to a negative color conductivity in 
the boost-invariant case. 



II. THE MODEL FORMULATION 

As was mentioned in Introduction, the classical Yang- 
Mills theory in space-time with pseudo-cylindrical metric 
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(r and 77 are proper time and space-time rapidity, re- 
spectively) has been abelianized since r rs 3/A s , where 
A s « 2 GeV It means that we actually come to the 
Maxwell theory with 4-potential (hereafter, we ne- 
glect the normalization constant l/\fW Cl where N c is the 
number of colors) . 

The free-field theory in mid-rapidity region in the case 
of central collisions, when the potentials are parametrized 
as A T = (CGC-like gauge fixing), A n = $(r, ry), 
A rT = 0, A v = $(r, rx), has been already examined 
(see Ref. (23[) in order to describe the space-time evolu- 
tion of the field flow (collective velocity) at pre-thermal 
stage of collisions. It turns out that the results obtained 
are qualitatively the same like in the case of non-Abelian 
model from Ref. 
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Here we generalize the abelianized Glasma equations 
by inclusion of sources in the right-hand side: 

d*<f> - -d T $ - d? T <Z> - — d rT $> = J„(r, r T ), (3) 
r tt 
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+ -d T V - cH * + — d rT V = J v {t, rr). (4) 

T r T 

Note that J T = 0, J rT = and the current conservation 
is satisfied automatically. 

In the context of A+A collisions, the presence of the 
sources corresponds to the existence of the medium. Ac- 
counting for hard partonic component (viewed as particle 
subsystem), we aim to investigate the field and particle 
dynamics. 

The components of the current in Minkowskian space- 
time are 



J" =9 M/+-/-) 



d 3 p 



(5) 



where p° = |p| (the case of massless partons), f± are 
distribution functions of (scalar) partons. The distribu- 
tion of partons is supposed to be anisotropic in the mo- 
mentum space and inhomogeneous in configuration one. 
Space-time development of functions f± is determined by 
Vlasov equations which we will formulate below. Note 
that "—g" corresponds to the charge of electron in the 
context of electrodynamics. 

A toy field model with non-trivial right-hand side has 
been already investigated in Ref. [26j . 

It is useful to parametrize momenta as (j/ 1 ) = 
(pr cosh y, pt cos <j>, pt sin <j>, pt sinh y) , where y is mo- 
mentum rapidity. 

In the terms of our variables, one has: 
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p T cosh8(f + - f_)dp T dydcj>, 
-T9 J Pt sinh0(/ + - f-)dp T dyd^>, 
z -9 J Pt C0S C(/+ - f-)dp T dyd<p, 
-r T g / p T sm£(f + - f-)dp T dyd(j), 
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where 8 — y — T), £ = <fi — ip. 

Taking conditions J T — 0, J rT — into account, the 
difference /+ — /_ should be odd function of 8 and £ 
during evolution. 

The evolution of f± is generated by Vlasov equations: 

(L±gF)f±=0, (10) 

where L = p^, F = pPF^, F„ v = d„A v - d v A^. 

Since the sources (partons) are randomly distributed 
at the initial moment to, we put /? = /_ = / , where 
/° is defined as (dNh/d 3 xd 3 p)\ To , and Lf° = in our 
investigations. It means that the system is neutral and 
the currents are absent at To. 

Using the curved coordinates, we obtain 



L = pt [ cosh 8d T 
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-(sinh 8 singdy + cosh 8 cos £8$) H — d$. (12) 



We can see that the operator of Lorentz force F is simply 
the operator of rotation in momentum space and, there- 
fore, conserves the absolute value of transverse momen- 
tum pt- It is expected that such structure of the Lorentz 
force should lead to the momentum transmission between 
different directions and, consequently, to the instabilities 
in this system. 



III. SOLUTION OF EQUATIONS 

In this Section, we concentrate on finding the solu- 
tion of the set of the coupled Maxwell- Vlasov equations. 
Since the Glasma field is essential at early stage of nuclear 
collisions (in contrast with the hard partons or quarks), 
the study of the field dynamics actually dominates. By 
this way, it is necessary to express the particle currents 
through the fields. In general, the Vlasov equations are 
complicated. For this reason, we are forced to use a 
method for approximate solving this set. 

Fluctuation of distribution function, which arises dur- 
ing fairly small time interval At, can be found in the 
linear approximation in g: 



f± = f°TgSf. 



(13) 



In this approximation, the space-time evolution of cor- 
rection 5f, determining difference / + — /_ = —2g5f and 
the current components, results from the following equa- 
tion: 



LSf = Ff 



(14) 
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Note that this approximation does not permit us to 
investigate an isotropization of the particle (hard par- 
ton) kinematic part of the energy-momentum tensor, pro- 
portional to the sum / + + /_ . It is expected that such 
isotropization effect can be observed if the correction of 
the order g 1 is included. Nevertheless, the approximation 
under consideration allows one to study the instabilities 
in the system. 

If r — ¥ To, there exists an approximate solution, which 
is short-living in time and localized in space, 



Sf 



Pt cosh 9 



Ff. 
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It is easy to verify that the action of evolution operator 
L on this expression gives us 



Pt 
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where 



W{t) = 



tanh 6 , r 

d T + (a 



coshf? TT 



•„ -r tanh 6 
tt cosh 9 v 



(Ff). (17) 



Thus, if t — > tq, the last term in right-hand side of 
(fTB")) vanishes. Also note that this is the simplified proof 
of Eq. (|15p . In order to understand this approximation 
in details, see Appendix A. 

Often the model initial boost-invariant distributions 
in central heavy-ion collisions take the form / = 
f°(pT,6) = f°{PT,—9) (note that for the sake of cor- 
rectness f° has to be also a function of rr). In this case, 
we obtain 



J n = (J- n {T)d T <S>, 



J v = (Tp(r)9 T *, 



(18) 



(dNh/d 3 x)\ To (rT) in transverse plane. Taking into ac- 
count that the initial distribution is even function of 8, 
one gets 
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where A is a positive constant arising after integration 
over momentum variables. 

Thus, <r„ < 0, while > 0. It means that the color 
negative conductivity takes place in the transverse plane. 

At this stage the natural question arises: how the 
negative conductivity does look in the laboratory ref- 
erence frame. Eqs. (TT81 are actually the Ohm law, 



where E„ 



± rr/ 



are the (chromo)electric field strengths. Introducing 
Ei — F t i in the Minkowskian space-time, we find that 
E n = tE z , E v = tt cosh?y(— sinipE x + cos tpE y ) + 
tt sinh rj(— sin <*pF zx + cos (pF zy ). In these terms the cur- 
rent components are 

Jt = — smhrjo^Ez, J z = cosh rja^Ez, (22) 
_ siny _ cos ip 

If rj = and ip — 0, one has that J t = J x = 0, J y — 
a v E yi J z — <7 r] E z . Thus the color negative conductivity 
takes place under some conditions (related with the value 
of angles) in the laboratory reference frame. This effect 
is associated with filamentation in the plasma [3J|. 

Since it is hard to find the general solution of field 
equations for arbitrary distribution (dNh/d x)\ To , we try 
to study the particular case, when (diV^,/c? 3 a;)| T0 = const. 
This assumption simplifies the problem significantly. 

When cr is a constant, the spatial dependence of 
the field potentials is immediately derived by using the 
Bessel-Fourier transform: 



where 



(t v (t) = 2(r - T )cr , ov(t) = -( r - T )a (19) 

are conductivities. 

The common multiplier dependent on the initial dis- 
tribution of partons is 

oo oo 

(7 = -27r.g 2 J dp T J dyd y fpT tanh 6. (20) 



We can immediately see that a v (r) and <j v {t) have 
the different signs. It says about the presence of negative 
color conductivity driving to instability in the system. 
The mechanism of this instability looks simple: we deal 
with situation when the particles (partons) give the en- 
ergy to the field. 

Now let us analyze the properties of o~$. Firstly, we 
assume that the initial distribution / is the product 
of the function of (pt,0) and the spatial distribution 



$(r,r T ) = / <f>o(kr) gT,(T,k T )MkTrT)dk T , (24) 



OO 

*(r,r T ) = r T J * (/c T ) 5v (r, k T ) Ji(fc T r T )cWc T ,(25) 



where the initial conditions resulting from CGC concept 
are applied: 



9v\ra = 0, 



= fc 7 



9 v \r = I, Td T g v \ ro = 0. 



(26) 
(27) 



In principal, functions ^(r, kr), ^(t, fcr) can be ex- 
pressed for arbitrary tq =/= in terms of Heun functions. 
However, these expressions are complicated for heuristic 
analysis of our model and its applications. 

For this reason, we write down the functions g v v at 
to -> 0: 
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where M(a, b; z) is the Whittaker function. 

It is easy to verify that the occurrence of negative 
color conductivity a v leads to a growth of the some com- 
ponents of magnetic and electric fields in comparison 
with the case of the theory without partonic medium. 
It is important that the Abelian magnetic field exhibits 
a growth, draining some energy from the particle reser- 
voir. The instabilities are related with the presence of 
exponents in functions gr),tpi Whittaker functions change 
actually a phase of oscillations only in comparison with 
the free theory. 

If <7o — ► and tq — > 0, we come to the well-known 
expressions (see, for example, Ref. [27j and references 
therein) : 

g v (r, k T ) = TJi(krT), g v {r, k T ) = J (/c T r), (29) 

where J n (z) is the Bessel function. 

These expressions correspond to the perturbative (low- 
est order in the source charge densities) solution. 

Now it is necessary to determine the functions <I > o(^t) 
and ^o(^t)' They originate from the initial conditions 
for the field equations. 

Note that ^o(^t) and ^(^t) are fluctuating quanti- 
ties within the CGC concept, and the pair correlator of 
the (Yang-Mills) potentials is observable only. However 
the field potentials in our approach are not stochastic 
quantities in the contrast with CGC ideology because 
we goal to constitute the initial conditions on the base 
of the statistically averaged components of the energy- 
momentum tensor accounting for the spatial inhomo- 
geneity. 

In order to derive ^o(kr) and $o(&t), let us use the en- 
ergy density distribution and the requirement of the ab- 
sence of field flow at the initial moment. In mid-rapidity 
region (77 = 0) and ip = (note that transverse directions 
are equal in the system with cylindrical symmetry) , when 



Ttt — T Tr , T tx 



T TX , we have that 



Ti 



tt\r 



T, 



tx\T 




(30) 
(31) 



where s{tt) is assumed to be the known function from 
numerical calculations or physical point of view. 

Our trick consists in division of the energy density be- 
tween different field components: 



*|r = Var T f(r T ), - J — 

T 



= Vl - af(r T ), 



where /(fr) = y^i^Tr) and a is a some separation con- 
stant (in general, a should be function of rr). Since the 
potentials Vt, $ are real, one has that < a < 1. 

In principal, a is arbitrary constant. Practically, it 
turns out that awl/2 (it follows from comparison of 
electric and magnetic strengths within the numerical ap- 
proach). Note that the observables of the source-free 
theory are independent on a. 

Thus, one finds 



*o(M = V^f(kr), $o(fcr) = VT^ af(k T ), (33) 
here 

00 

f{kr) = J f{r T )Jo{k T rT)rTdr T . (34) 


These expressions finally determines the fields in our 
model. 



IV. APPLICATIONS 

In the previous Sections we have formulated the model 
of Glasma in hard partonic medium. Since the classi- 
cal field modes are usually interpreted as soft partons, 
their spatial dependence at early stage of A+A collision 
may be done within the framework of Glauber model. 
However, the explanations of experimental data can be 
efficiently done with application of another distributions 
too. As it was demonstrated in Ref. [24], the Gaussian 
distribution of soft partons leads to the adequate pion 
spectra produced after collision at RHIC. To formulate 
the field initial conditions, here we would like to choose 
the same approximation for the energy density at the 
initial moment, 



s{tt) — E exp 

where E = 45 GeV/fm 3 , R 
Then we find that 



2R 2 1 



3.768 fm. 



f{k T ) = 2 3/2 VER 2 exp (-fcf-i? 2 ). 



(35) 



(36) 



(32) 



The boost-invariant distribution function / (defining 
conductivities) is completely arbitrary at this point, so 
in order to proceed one needs to assume a specific form 
for it. In what follows we will require that / is obtained 
from isotropic function, 



(37) 



N exp 

' Ph 
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by the replacement y — > 6 in p° = px cosh y and by the 
rescaling of one dimension in momentum space, 



/° = iV(C)exp 



Pt 
Ph 



cosh 2 + ( sinh 2 I 



(38) 



where takes the role of saturation moment, £ > — 1 is a 
parameter reflecting the strength of the partonic medium 
anisotropy and N(£) is a normalization constant. Note 
that C > corresponds to a contraction of the distribu- 
tion in the z-direction whereas — 1 < £ < corresponds 
to a stretching of the distribution in the z-direction. 

Constant iV(£) is simply determined by requiring the 
number density to be the same both for isotropic and 
anisotropic systems and can be evaluated (by integration 
over momentum variables) to be 



JV(0 = No^/T+(. 



(39) 



Integrating over momentum, the multiplier defining 
the conductivities is 



a = Airg 2 N piC{Q, 



where 



C(C)= ga + C) 372 ^ 



4 



< 



(40) 



(41) 



The coefficient C(£) in the region £ e (— l,oo) is de- 
termined by hypergeometric function F and is such that 
C(— 1) = (the case of source- free theory), C(0) = 2/3 
(for isotropic medium), C(oo) = tt/2. 

Fig. 1 shows how the exponentially growing energy 
transferred from hard to soft partons is distributed 
among magnetic and electric fields at ao ^ in com- 
parison with the case of free field theory. The dominant 
contribution is still in longitudinal electric field (in accor- 
dance with CGC-like initial conditions). Nevertheless, we 
see that the transverse magnetic field demonstrates un- 
stable behavior too while this effect is absent in the free 
theory. Since the particle subsystem gives the energy to 
the field, the total field energy density (as the sum of 
components) tends to grow with time. 

Note that a similar model with expanding Abelian field 
coupled to the hard partons, when the strict boost invari- 
ance of fields is relaxed, has been already developed in 
Ref. [28| in the context of the quark-gluon plasma. 




FIG. 1: Time evolution of the field energy density split into 
longitudinal and transverse electric (E) and magnetic (B) 
components at rr = 0.1 [fin], r\ = 0, tp = 0. The top panel 
corresponds to the free theory with <to = 0. The bottom panel 
demonstrates the growth of El and Bt at ctq = 0.25 [fm -2 ]. 



contrast with the usual approach based on the energy- 
momentum tensor. As the result, the instabilities in 
the system under consideration lead to a conclusion of 
the presence of negative color conductivity in A+A col- 
lisions. The sign of conductivity in transverse plane de- 
pends on the angle what says about filamentation inher- 
ent to Weibel instabilities in plasma. 

Unfortunately, the approximate solution derived to the 
transport equations does not permit us to achieve the 
isotropization here. This problem should be investigated 
in details and will be published elsewhere. 



V. CONCLUSIONS 



Generalizing the space-time evolution of the expand- 
ing Glasma with CGC-like initial conditions by inclusion 
of small density of the hard partons anisotropically dis- 
tributed in momentum space, we observe the instabilities 
(due to transferring energy from hard to soft partons) in 
the case of the abelianized boost-invariant model. 

Here we propose to measure an anisotropy by means of 
transport coefficients like the conductivity tensor in the 
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APPENDIX A. PROPAGATOR OF TRANSPORT 
EQUATION 

Here, we would like to discuss in details the approxi- 
mation which we have applied previously. To find 8f, we 
have to determine operator L^ 1 , inverse to the first-order 
evolution operator L. The inversion procedure of the evo- 
lution operator of transport equation was elaborated by 
Landau and, generally speaking, results in emergence of 
the Landau damping in plasma. 

Let G is the solution of the following equation: 



LG{r,e,v T \r',9'y T ) 



5{t - t') 



5{e-e')5\v T -v' T ), 

(42) 



where right hand side is the 4-dimensional 5-function 
with respect to the pseudo-cylindric measure rdrdOd ! r-r- 



Without loosing generality, we limit ourselves by the 
case, when r — > 0, and by using the transverse coordi- 
nates i"t instead of cylindrical (rr, f)- 



Using the formulas from Appendix B, G is represented 



as 



G(T,d,r T \T',9',r' T ) = lim 



+o (27r)Vr 



piuj[T cosh (£— 0)— t' cosh (£— B')\— ik T (r T — r' T ) 

k T VT — co cosh £ — ie 



uoduod^d 2 fcy, 



(43) 



r 



where v T = Pt/Pt and the Landau damping is already 
taken into account due to auxiliary formula: 



lim — - — = V— + i-k8{x). 

e^+a x — ie x 



(44) 



Discarding the spatial dispersion, we have to assume 
that kx <C cjcosh£. It leads to simplification: 

G(T,6,r T \T',6',r' T ) w — e(rcosh6» - r'cosh^') 
Pt 

x5(t smh6 - t' smh6')5 2 (r T - r' T ). (45) 

where 6 is Heaviside function defined as 9(a; < 0) = 0, 
Q(x = 0) = 1/2, Q(x > 0) = 1. 

This approximation says that the system is homoge- 
neous at significantly large times and the fluctuations 
are localized in space. 

Furthermore, let the Bjorken scaling flow, when 9 ss 0, 
take place. It means that 9 and 0' should be equal and 
gives us that 

G(r, 9, r T \r', 9\ v' T ) « ^^5(9 9>)5\v T - v' T ). 

(46) 

More precisely, it can be derived from condition, 
rsinh# = const, resulting in 



sinh 9dr + t cosh d9 = 0, 



(47) 



where dr = t - t' , dO = 6 - 6 1 . 

Assuming that r is small and the expression under in- 
tegration is not essentially changed at this time range, 
we can do the following replacement: 



oo 



Q{t-t')F{t') -> tF{t). 



(48) 



Then, one obtains that 

J G(t, 9, r T |r', 9\ t' t )F(t', 9', v' t )t' dr' d9' d 2 r' T 



Pt cosh 9 



F(t,9,t t ). 



(49) 



This formula determines the solution of inhomogeneous 
transport equation with a source F. 



APPENDIX B. INTEGRAL TRANSFORMATION 

The Fourier transformation reads 

oo oo 

/(«,*)= / |^ / fMeW-rt-^'-'Upda. 

— oc — oo 

Let us introduce new variables: 

t = rcosh9 7 z = rsinh#, 
p = pcoship, s — psmhtp, 
fi = A cosh (j), fi' = A sinh (f>. 

If f(t, z) — F(t, 9), one gets the following transforma- 
tion rule: 



F(r,9) 



OO 

•*»-/ 

— OO 

oo 



F(p,V)e- 



-ipX cosh (4>—ip) 



F(A,0)e 



ir\ cosh [<j>— 9) 



pdpdtp, 



Xd\d(j>. 
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For example, we find that where At = r — tq, A(9 = 6 — 8q. 

oc 

t (2tt)^ J 

— oo 
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